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■ Abstract 

O ■ 

' We examine the deformation quantization of a single particle moving in 

I ' one dimension (i) in the presence of an infinite potential wall, (ii) confined 

, by an infinite square well, and (iii) bound by a delta function potential en- 

f*^ ' ergy. In deformation quantization, considered as an autonomous formula- 

tion of quantum mechanics, the Wigner function of stationary states must 
be found by solving the so-called *-genvalue ( "stargenvalue" ) equation for 
the Hamiltonian. For the cases considered here, this pseudo- differential 
equation is difficult to solve directly, without an ad hoc modification of 
the potential. Here we treat the infinite wall as the limit of a solvable 
r-{ ' exponential potential. Before the limit is taken, the corresponding *- 

genvalue equation involves the Wigner function at momenta translated 
IL" ' by imaginary amounts. We show that it can be converted to a partial 

. | differential equation, however, with a well-defined limit. We demonstrate 

■ that the Wigner functions calculated from the standard Schrodinger wave 
' functions satisfy the resulting new equation. Finally, we show how our 

results may be adapted to allow for the presence of another, non-singular 
part in the potential. 
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1 Introduction 



Deformation quantization 1 is touted as a completely autonomous method of 
doing quantum mechanics, that is especially suited to the study of the classical 
limit. 

Some quantum systems that are simple to treat in the Schrodinger formu- 
lation, however, are difficult to analyze in deformation quantization. Even the 
treatment of a free particle is not completely straightforward. 2 Other examples 
are a particle moving in one dimension in the presence of an infinite wall, and 
a particle confined by an infinite square well @|. Surely, if deformation quanti- 
zation is to take its rightful place as one of the possible ways of doing quantum 
mechanics, such simple systems must be treatable in it. 3 

For stationary, pure states, the *-genvalue ( "stargenvalue" ) equation must 
be solved to find the Wigner function. We will restrict to such considerations 
in this work. For the infinite wall(s) case, simply imposing the usual boundary 
conditions on its solutions does not lead to the expected Wigner functions. In 
0] it was found that if the infinite wall and infinite square well potentials were 
modified by an additional "boundary potential" in an ad hoc way, then the 
expected Weyl transform of the density operator was a solution. The authors 
of @] did also show that the added potential terms were consistent with the 
Schrodinger treatment of those systems. They did not, however, derive the 
terms from first principles. 

Furthermore, the normal intimate relation of deformation quantization with 
classical mechanics is altered by their modification. In the deformation quan- 
tization of standard systems, quantum mechanics is treated as a deformation 
of classical mechanics, with deformation parameter h. The Hamiltonian is not 
modified by quantum corrections, and the deformation is encoded entirely in 
the ^-product. On the other hand, the "boundary term" added to the potential 
in 0] is of order h 2 , taking the form —^-8'{x — a), where x — a is the position 
of an infinite wall. 

Here we study the one-dimensional infinite wall by treating it as the limit of 
a solvable exponential potential, limQ^oo e 2ax . Our original hope was to derive 
the prescription of jl] from first principles. We follow instead a different path, 

1 Deformation quantization is also known as the Weyl-Wigner(-Moyal) formalism, phase- 
space quantization, and by other names as well. For elementary introductions, see [Q. More 
advanced reviews are listed in [2]. The modern interpretation of the formalism as a deformation 
of classical mechanics was first reported, and developed extensively, in 

2 See our Appendix. 

3 MW thanks Brian Wynder for emphasizing this point to him. 
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however. 

Before the limit of the exponential potential is taken, the *-genvalue equa- 
tion involves the Wigner function evaluated at momenta shifted by imaginary 
amounts. Finding the limit of this equation directly is problematic. We show, 
however, that it can be converted to a partial differential equation, with a well- 
behaved limit. The resulting new equation can be treated in the standard way: 
no modification of the potential is required, and the normal Dirichlet boundary 
conditions can be used. We demonstrate that the Wigner functions calculated 
from the standard Schrodinger wave functions satisfy the new equation. Finally, 
we adapt our results to allow for the presence of another, non-singular part in 
the potential. 

We should point out here that even the Schrodinger (canonical, operator) 
quantization of infinite walls has some subtleties (see |3] and references therein) . 
In their operator quantization of such one-dimensional potentials, the authors 
of |S] adopt an approach similar to ours: they resolve apparent paradoxes by 
"acknowledging the existence of the rest of the real line." Their treatment has 
no need for additional boundary potential terms like those introduced in 0]. 

A quick review of deformation quantization is given in section 2. Section 3 
is the meat of the paper, containing the treatments of various cases of infinite 
walls and wells described as limits of exponentials. In section 4, the results 
are adapted to include an additional, regular potential, besides those described 
in section 3. Section 5 is our conclusion, and the Appendix treats the pure 
deformation quantization of a free particle in one dimension. 

2 Deformation quantization 

This section is a quick review of deformation quantization ^ [2], setting our 
notation and providing the results we will need. We will restrict attention to 
pure states, described by a state \ip), or a density operator p — \ip)(tp\. 

In deformation quantization, observables are not represented by operators, 
but rather by functions on phase space. They are multiplied using a pseudo- 
differential ★-product ("star product") that is associative but non-commutative. 
Included is the Wigner function describing the state of the system, the central 
object in deformation quantization. It obeys an evolution equation involving a *- 
commutator. In the h — > limit, *-commutators reduce to Poisson brackets, and 
the equations of motion of classical mechanics are recovered. In that sense, this 
autonomous method describes quantum mechanics as a deformation of classical 
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mechanics 3 . 

Deformation quantization can be understood as a transform of the standard 
way of doing quantum mechanics using the density operator (matrix). Let us 
restrict attention (throughout) to the case of a single particle moving on the 
x-axis, so that phase space has coordinates (x,p). 

In canonical quantization, the phase space coordinates x and p are pro- 
moted to the operators x and p, obeying the Heisenberg commutation relation 
[x,p] = ih. Observables like x 2 p are promoted according to an operator order- 
ing scheme. 4 Choosing the Weyl ordering, our example becomes 

W~ 1 (x 2 p) = ]-(x 2 p + xpx+px 2 ) . (1) 
For a function f(x,p) on phase space, this generalizes to 

W" 1 / - f(d a ,d b )e a£+b %^ , (2) 
where d a := J^, etc. In particular, it follows that 

This gives rise to another expression 

W ~ lf = (2^)2 / drdadxdp f(x, P y T( -P- p ') +l < T( - & -^ , (4) 

where (01 has been used in the usual Fourier formula for /. 
The crucial property of this ordering is 

(W- 1 /) (W-\ g ) = W-\f*g) , (5) 

proved by Gronewold. The Gronewold-Moyal ^-product (pronounced star-product) 
takes the form 

f(x,p)*g(x,p) := f(x,p) expjy ( d x d P - d P d x )} g{x,p) . (6) 

Here the arrows indicate the directions in which the derivatives act, and they 
act only on / and g, and not to the left or right of them. That is, eqn. |JBJ 
stands for 

f(x,p)*g(x,p) = e* f( x ',p')g(x,p) | a , =X| p>=p . (7) 



4 Deformation quantization may be carried out in ways that correspond to many different 
operator orderings. Here, we will restrict to the Weyl ordering W _1 , resulting in the famous 
Gronewold-Moyal *-product. The use of the inverse notation (W -1 ) is for later convenience. 
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Eqn. JSJ) shows that the ^-product is an associative, non-commutative prod- 
uct of phase space functions that mimics the product of operators. In deforma- 
tion quantization, the operator product is replaced by this ^-product, so that 
operators can be represented by ordinary functions (and distributions) on phase 
space. 

The exponent of the ^-product JSJ) indicates the most important property of 
deformation quantization: its intimate relation to classical physics. In classical 
mechanics, it is the Poisson bracket of functions on phase space, 

that enters the dynamical equations. In the operator formulation of quantum 
mechanics, it is the commutator [/,<?] of operator observables / and g that is 
important. In deformation quantization, the ★-commutator 

[f,gl := f*g-g*f (9) 

of functions / and g takes its place. The equation 

limi[/,g]* = {/,<?} (10) 
h^o in 

encodes the correspondence relation between classical and quantum mechanics 
in deformation quantization. 

The inverse W of the Weyl map W _1 is known as the Weyl transform. One 
formula for it, 

w f = _L f d^dip Tr(f e *K0-p)M-(*-*)*>]/ft) ( ( U ) 
2i:h J 

treats momenta and coordinates in symmetric fashion. A useful variant can be 
derived from it: 

Wf = % j dye-*P«(x+ ^|/>- ^) , (12) 

as well as a similar formula involving momentum eigenstates. The Weyl trans- 
form Wf of the operator / is also known as its Weyl symbol. 

The object 

A(x,p;x,p) := e *[tf-P)C+(«-*)<p]/1 , (13) 

is prominent in the fundamental equation (jl 1|) . It is sometimes called the quan- 
tizer, and its importance was stressed in (i . 

The Weyl transform obeys 

W(f)*W(g) = W(fg) , (14) 
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an inverse analog of JSJ. This makes it possible to work exclusively with phase 
space functions, as long as they are multiplied with the ^-product. Deformation 
quantization is the realization of this idea - it is the Weyl transform of quantum 
mechanics done with the density operator. 

The Weyl transform of the density matrix p = \ip) {ip\, 

W{ P ) = h fdye- i ^r(x-^-)Hx + ^-), (15) 

is the central object in deformation quantization; it describes the quantum state 
of the system. After normalization, it is known as the Wigner function: 

P : = ■ (16) 

Combining the last two equations gives 

*#] := ^ J dye-^r(x- T f)iix + ^-) . (17) 

We will use the notation p[ijj\ to emphasize that this is the Wigner function 
calculated from a known Schrodinger wave function rp. 

The Wigner function evolves according to 

ih^ = [H,p\*. (18) 

For stationary states, ^| = 0, so that 

[B,p\* = . (19) 

For p describing an energy eigenstate with eigenvalue E, 

H*p = p-kH = Ep . (20) 

These simplified dynamical equations allow us to solve for the Wigner function 
of the stationary states of the system. 

By <|15|) . the Wigner function can always be calculated if the Schrodinger 
wave function is known. We are interested here, however, in considering defor- 
mation quantization as an autonomous method of doing quantum mechanics. 
That means we want to find the Wigner functions by solving either (I18|) or (|20() , 
and we will restrict to the latter in this paper. Of course, l|15f) can still be used 
as a very useful check of our results. 

Solving l|20[) can lead to more general solutions than those of physical inter- 
est. Imposing the constraints 

p-kp = j^p , p* = p , (21) 



G 



is important in finding the correct solutions [3]. Clearly, the first relation is just 
the Weyl transform of the usual projection condition on the density operator. 
More generally, the requirement is 

Pi *Pj = v Pj , (22) 
for Wigner functions describing a discrete set of states, and 

Pa*Pp = j^S(a~f3)pp (23) 

for states with a continuous label. 

All observable probabilities can be calculated using the Wigner function. 
First of all, by JTSJl, the probability densities are 

mx)\ 2 = fdpp(x,p), m?)\ 2 = I dx p{x lP ) . (24) 



Clearly then, the Wigner function is normalized and real: J dxdpp — 1 and 
p* = p. The expectation value of an operator 

</> = Jdxdpp*f , (25) 

where / := W(/). Roughly, one can think of the integral over phase space as 
the analog of the trace, and as discussed above, the star product takes the place 
of the operator product. The important cyclic property of a trace is encoded in 

Jdxdpf-kg — Jdxdpg*f = Jdxdpfg. (26) 



3 Free particle except for infinite walls/wells 

In this section, we will treat a particle moving freely on the £-axis, except for 
the presence of one or more infinite potential walls or wells. More precisely, the 
infinite barrier, the infinite square well, and the delta-function well will all be 
studied as limits of potentials built from exponentials. 

We are interested in the "pure" deformation quantization of these systems. 
That is, we would like to derive their Wigner functions using the equations of 
deformation quantization only, without reference to the well known Schrodinger 
wave functions, for example. Here we restrict consideration to stationary states. 
That means, therefore, that we must examine the *-genvalue equation l|2(JII for 
the corresponding Hamiltonians. 
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As mentioned in the introduction, these problems are not straightforward, 
as we now explain 0]. Consider equation (|2U[1 more carefully, for the example 
of an infinite wall, i.e. for the potential energy 

V(x) = { °' X< °> (27) 
w oo, x > 0. y ' 

Following the Schrodinger treatment of this system, one would restrict to x < 0, 
and impose the boundary condition p(0,p) = 0. For x < 0, the *-genvalue 
equation is that of a free particle, with real and imaginary parts given in the 
Appendix as equations (|53|l and respectively. 5 But the imaginary part, 

p d x p = 0, does not lead to sensible results for this potential (see the Appendix). 

To look for some guidance, we can study the expected solution, the Weyl 
transform of the known density matrix. The Schrodinger wave function is 

tP(x) = 6{-x) \e lV ^ x - e~^ x ] . (28) 

Using (|17f) . the corresponding Wigner function is found to be 

p[0 = 9(-x)p(x,p) , (29) 

with 

s _ 2 S in[2x(p+y^E)] 2 sm[2x(p-VE)] 

+ 2cos(2xV^) 2si "( 2 *P) . (30) 

This p(x,p) does not satisfy the imaginary part l|84|) of the *-genvalue equation 
H -k p = Ep. 

Our goal in this article is to find a dynamical equation replacing the *- 
genvalue equation that can be solved to find the correct Wigner functions for 
these problems. 

3.1 Infinite wall 

The Liouvillc Hamiltonian is 

H a = P 2 + e 2ax . (31) 

Constants like the mass have been set so as to simplify considerations (2m = 1, 
e.g.). A pure deformation quantization of this system has already been carried 

5 In those equations and henceforth, we set h = 1. 
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out in 7 . Taking the a — > oo limit of V(x) = e 2ax yields an infinite wall with 
V = for x < 0, and = oo for x > 0. 

Therefore, the limit of the solution found in |7| yields the correct Wigner 
function. That is not what we wish to do, however. We hope instead to take 
the limit of the *-genvalue equation, in order to find an equation that can be 
solved directly, leading to the physical Wigner function. 

The *-genvalue equation is easily found to be 



H a * p(x,p) 



p{x,p) = Ep{x lP ) . (32) 



It separates into S (imaginary) part 

[-pd x + e 2ax sm(ad p )]p(x,p) = 

and Sft (real) part 

1 



p* - E - j dl + e 2ax cos(a d p ) 



Formally, these equations can be rewritten as 



p(x,p) = . 



-2ax 



d x p{x,p) 



2p 



[p(x, p + ia) — p(x, p — ia)} 



(33) 



(34) 



(35) 



and 



p 2 -E-]d 2 )p(x,p) 







(36) 



Using (|35f) to find <9^p(#,p) ; and substituting this into l |36f k we arrive at an 
equation without derivatives, 



p(x ,p-\-2ia) — p(x ,p) . p(x,p— 2m) — p(x,p) 
p-\-iai p—ia 



= ^-E)p{x,p) + \{f^) 

Ip( x ,P + ia ) " p(%,P~ ia)] 
+ ^2~ [p( x ,P + i a ) + p{%iP- ioi)] 



(37) 



This can be considered a difference equation in the momentum variable. Only 
imaginary shifts of the momentum arguments are involved; that is, besides 
p(x,p), the result involves the four quantities 



p(x,p±ia) , p(x,p±2ia) 



(38) 
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The a — > oo limit of (|37[) is problematic. 

We can trade the four quantities of however, for the derivatives 

a%p{x,p) , n=l,2,3,4. (39) 

The resulting differential equation will have a well-defined limit as a — > oo, the 
result we are seeking. 6 

Four equations relating the "variables" of (|38|l to those of l|39[) are required. 
Two are already provided: i|35|) and l|36|l . The two additional equations can be 
derived by taking derivatives of 136(1 : 

= dlp(x,p)-4(p 2 -E)d xP (x,p) 

+ Aae 2ax [p(x,p + ia) + p{x,p — ia)] 
- 2e 2ax [d x p(x, p + ia) + d x p(x, p - ia)] , (40) 

and 

= d i p(x,p)-4(p 2 ~E)d 2 xP (x,p) 



8a e [p(x,p + ia) + p(x,p - ia) 



4aie 4ax 



p(x,p + 2ia) — p{x,p) p(x,p — 2ia) — p(x,p) 



p + ia p — ia 

2e 2ax | [4(p + ia) 2 - E]p{x,p + ta) 

+ 2e 2ax [p(x,p + 2ia) + p(x,p)] 
+ [4(p — ia) 2 — _E] p(x,p — ia) 

+ 2e 2ax lp(x,p-2ia)+p(x : p)]} . (41) 



With the help of symbolic computation, we find a simple, differential equa- 
tion results: 

= ±- Q dip{x,p)+ l -(p 2 + E)d 2 xP {x 1 p) 

+ (p 4 - 2Ep + E 2 ) p(x,p) - e iax p(x,p) . (42) 

Taking the limit a — > oo, we find the new equation 

^ dtp{x 7 p) + l - (p 2 + E) d 2 x p{x,p) + (p 4 - 2Ep + E 2 )p{x,p) = , (43) 



6 Holder's theorem (see e.g.) states that no solution of the simple difference equation 
y(x + 1) — y(x) = 1/x satisfies any algebraic differential equation. Replacing a difference 
equation with a differential equation may, therefore, omit interesting solutions. In our case, 
however, we will see that the physical Wigner function is a solution of our resulting differential 
equation. 
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valid for x < 0. 

It is simple to verify that the Wigner function p{x,p) of (|30fl satisfies the 
new equation (|43() . for x < 0. 

It is also interesting to note that eqn. I|43|) can be rewritten as 

(p 2 * p*p 2 - E 2 p) - 2E$.(p 2 * p - Ep) = 0. (44) 

The new equation is a linear combination of ones that follow from the -Ar-genvalue 
equation for a free particle. 



3.2 Infinite square well 

We will now consider the case of two infinite walls, as limits of two exponential 
potentials. More precisely, we will be taking the a — > oo limit of the *-genvalue 
equation following from the sinh-Gordon Hamiltonian 



H,, 



p 2 + e-M^+i) + e M*-i) 



(45) 



The method and the result for the two-wall potential will be remarkably similar 
to those for the single-wall case. 

The *-genvalue equation is 



It has 3 and 3? parts 



2a ( x + — d p 



h 2e~ 2a cosh 
= H a *p(x,p) = Ep(x,p) 



p(x,p) 



d x p{x,p) 



-2a 



[p(x,p + ia) — p(x,p — ia)] sinh(2o!a;) 



(46) 



(47) 



and 



p 2 -E--dljp(x,p) 
+ e 



[p(x,p + ia) + p{x,p — ia)] cosh(2ax) = 0. (48) 



Combining the previous two equations leads to 
= (p 2 -E)p(x lP ) 



4p 



■ cosh 2 (2aa;) 



p(x + 2ia) — p{x,p) p{x,p) — p(x,p — 2ia)) 



p + ia 



p — ia 



2p 



■ cosh(2ax) [p(x,p + ia) — p(x,p — ia)] 



-2a 



cosh(2aa;) [p(x,p + ia) + p(x,p — ia)] 



(49) 
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This equation involves no derivatives, but only the quantities 

p(x,p±ia), p(x,p±2ia) 



(50) 



as well as p(x,p). We again wish to eliminate these in favor of derivatives <9™p, 
n = 1, 2, 3, 4. To do so, we make use of the additional relations 

dlp{x,p) = 4(p 2 -E)d x p(x,p) 

+ 2ae~ 2a sinh(2aa;) [p(x,p + ia) + p{x,p — ia)] 
— 4ie~ 4Q cosh(2aa;) sinh(2o:a;) x 



p(x,p + 2ia) - p(x,p) p(x,p) - p(x,p-2ia) 



p + ia 



p — ia 



(51) 



and 



dtp(x,p) = 4(p 2 - E)d 2 p(x,p) 

+ I6a 2 e~ 2a cosh(2ax) [p{x,p + ia) + p(x,p — ia)] 
+ 16iae~ Aa sinh 2 (2aa;) x 

p{x ,p+2ia) — p(x.p) . p(x.p) — p{x,p— 2ia) 
p+ia p—ia 

+ e~ 2a cosh(2aa;) | [4(p + ia) 2 — E\ p{x, p + ia) 

— e~ 2a cosh(2aa;) [p(x,p + 2ia) + p(x,p)] 
+ [4(p — ia) 2 — E\ p(x,p — ia) 

— e~ 2a cosh(2as) [p(x,p) + p(x,p — 2ia)] 



(52) 



The resulting new equation is quite complicated, and so we refrain from 
writing it here. It does, however, have a well-defined limit a — > oo, for x S 

]-Mh 

^ dip{x, p)+ l - (p 2 + E) d 2 p(x,p) + (p 4 - 2Ep + E 2 )p(x,p) = . (53) 

This is identical to the analogous result l|43|l found for the one-wall case. 

As for the one- wall case, the Wigner functions calculated from the Schrodinger 
wave functions 



tp(x) = 9{-x + l)6(x + 1) cos(V-Ea;) , E = 



(54) 



satisfy this new equation, in the following sense. Using (|17f) . we find p[ip] 
6{-x + l)6(x + 1) p, where 

sin[(2p + mr)(l - \x\)] sin \{2p - rm) (1 - |a:|)] 



p(x,p) = 



2p + Tin 



2p — nn 
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or, equivalently, 



sin 



p(x,p) = 



2(p+s/E)(l-\x\ 



sin 



2(p-VE)(l-\x\) 



2{p + VE) 2(p - VE) 

cos(VEx) sin[2p(l- \x\)] 



(56) 



p satisfies the new equation l|53|). As in the one- wall case, the expression for the 
Wigner function valid where V — solves the derived equation. 

3.3 Delta-function potential well 

A ^-function potential well can be studied as the a — > oo limit of the Hamiltonian 

H a = p 2 ~ 2ae- 2a ^ . (57) 
Restricting to x > 0, the *-genvalue equation is 

H a *p(x,p) = 



p- l -d x ] -2ae- 2a ( x+ i d *) 



p{x,p) 



Ep(x,p) 



Its 3 and 9? parts are 



ict 



e ax d x p{x, p) = p + ia) - p(x, p - ia)] 



(58) 



(59) 



and 



e 2 ax [ 1 ?- E --diy{x 1 p) 

a \fl(x,p + ia) + p{x,p — ia))] = 



(60) 



Combining these last two relations yields an equation involving no deriva- 
tives, 

= 4(p 2 - E)p(x,p) - Aae~ 2ax [p(x + ia) + p(x,p-ia)] 

4a 2 e~ iax \ p(x,p + 2ia) — p{x,p) p{x,p — 2ia) — p{x,p) 



p p + ia p — ia 

2ia 2 e~ 2ax [p(x,p + ia) - p{x,p - ia)} , 



(61) 
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but instead depending on p(x,p) and 

p(x,p±ia), p(x,p ± 2ia) . 



(62) 



In the by-now familiar way, we trade the dependence on these for ^-derivatives, 
with the help of the additional equations 

dlp{x,p) = 4(p 2 - E)d xP {x,p) 

+ 8a 2 e~ 2ax [p(x,p + ia) + p(x,p-ia)] 



+ Aa 2 ie~ 4ax 



p(x, p + 2ia) - p(x, p) p{x,p) - p(x, p - 2ia) 



p + ia 



+ 



p — ia 



(63) 



and 



d 4 x p(x,p) = A(p 2 -E)d 2 x p(x lP ) 

- 16a 3 e~ 2ax [p(x,p + ia) + p(x,p-ia)] 



- I6a 3 e- 4ax i 



p{x,p + 2ia) — p{x,p) p{x,p) — p{x, p — 2ia) 



+ 



p — ia 



Aae 



p + ia 

[-4(p + ia) 2 + E] p(x,p + ia) 
Aae~ 2ax [p(x,p + 2ia) + p(x,p)]+ [-4(p - ia) 2 + E] p(x,p - ia) 
- 4ae- 2ax [p(x,p-2ia) + p(z,p)] J . 



(64) 



The resulting new differential equation 

^d 4 p(x,p) + ^ ^ ^ d 2 x p(x,p) 

+ [p 4 - 2Ep + E 2 + 4a 2 e- 4ax ] p{x,p) = 
has a well-defined limit. Taking a — > oo, gives 

-^d 4 p(x,p) + l -{p 2 +E)d 2 xP {x,p) + (p 4 -2Ep + E 2 )p(x,p) = 

since x > 0. The result is identical for x < 0, and to the result found above for 
both the one-wall and the infinite square well cases. 

The sole state bound by the delta- function well has wave function 



(65) 



(66) 



ip(x) = , E 



-1 . 



The corresponding Wigner function p[ip] is 

e -2x 

p(x,p) = 



cos(2xp) + i sin(2a;p) 



p 2 + 1 

It satisfies the new differential equation, where V = 0. 



(67) 



(68) 
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4 Wall with an additional potential 



So far, we have treated systems that are free except for infinite walls or wells. 
Describing such potentials as limits of exponential terms, we found in all cases 
that the same differential equation governs the Wigner functions where V = 0. 
It would be useful to generalize to systems involving additional, regular potential 
terms, in the presence of infinite walls/wells. 

As a first step, we will consider here the one-wall case, with an unspecified 
potential added. That is, we'll consider again the limit of the Liouville potential, 
with the additional, regular potential V(x) present. After a general treatment, 
we will consider the simple, special case V(x) = x 2 as a check. 

The Hamiltonian is 



H a = p 2 + e 2ax + V{x) 
yielding the *-genvalue equation 
H a *p(x,p) = Ep(x,p) 



(69) 



p- l -d x ) + e 2 < x+ i d >) \ p{ x ,p) + V(x) * p(x,p) . (70) 



Its 9 and 3? parts are 

[-pd x + e 2ax sm(ad p )] p(x,p) + 3 [V(x) * p(x,p)\ = 

and 



(71) 



p 2 - E - - O 2 + e 2ax cos(a d p ) 

Combining these gives 
= 4(p 2 -E)p(x,p) 



p(x,p) + K\y{x)*p(x,p)] - 0. (72) 



+ 



4p 



p(x + 2ia) — p(x,p) p(x — 2ia) — p(x,p) 



p + ia 



+ 



p — ia 



2e 



2ax 



[p(x,p + ia) - p(x,p-ia)] 

[p(x,p + ia) + p(x,p-ia)] 
1 



+ m[V(x)*p(x,p)] - -d x Z[V(x)*p(x,p)} 

p 

a 2ax ( ^[v(x)* p{x,p + ia)} 3 [V(x) * p{p - ia)} ' 



+ 



2p 



p + ia 



p — ia 



(73) 
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Following the procedure used above, we would like to eliminate the quantities 
p(x,p ± ia), p(x,p ± 2ia) in favor of derivatives of the Wigner function. As 
in the cases previously considered, two additional equations can be derived 
straightforwardly : 







d 3 x p(x,p) + 4(p 2 -E)d xP (x,p) 
+ 4ae 2ax [p(x,p + ia) + p(x,p-ia)] 



— ic 



p(x,p + 2ia) - p(x,p) pixjP) - p(x,p-2ia) 



p + ia 



p — ia. 



+ 2e 



2 (xx 



9 [V(x) * p(x,p + ia)] 3 W{ x ) * p( x ,P — ia 



+ 



p + ia p — ia 

+ 4d x ft[V{x)*p{x,p)} 



(74) 



and 

= 



d$p(x,p) + A(p 2 - E)dlp(x lP ) 

+ 8a 2 e 2ax [p(x,p + ia) + p(x,p-ia)] 

- 4iae 4ax 



p(x,p + 2ia) - p{x,p) pix^p) - p(x,p-2ia) 



p + ia 



p — ia 



+ 2e 2ax ^ [4(p + ia) 2 - E] p{x,p + ia) 
+ 2e 2ax [p(x,p + 2ia) + p{x,p)] 
+ [A(p-ia) 2 -E)] p(x,p-ia) + 2e 2ax [p(x,p - 2ia) + p(x,p)] J 

} 



+ 8ae 2 



2ax (%[V*p(x,p + ia)] $S[V*p(x,p-ia)] 



p + ia 



p — ia 



+ 8e Iax {&[V{x)*p{x,p + ia)] + [V(x) * p(x,p - ia)]} 

+ ±d 2 JK[V{x)*p{x,p)] . (75) 

This is not sufficient in this more general case, however. We now have 4 further 
quantities 



V{x) * p(x,p± ia) 



(76) 



to get rid of. Here (3J/3) indicates the real part or the imaginary part. 



With considerably more work, however, the required 4 equations can be 
derived. We find 



V(x)*[p 2 -E--d 2 p(x,p) 
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'2ax 



V{x) * p(x 1 p + ia) + p(x,p—ia) 



t/(x)*3? V(ar)*p(a;,p) 



(77) 



and 



(K/9f) 



V(a;) -kpd x p(x,p) 



— ($l/Z){V(x)*[p(x,p + ia) + p(x,p-ia)]} 
+ (31/3) jV(x)*3 V(s)*p(a;,p) } = 0. 



(78) 



By symbolic computation, the resulting differential equation can be derived. 
In the limit, it goes to 

^d$p(x,p) + <£±3.% p ( x , p ) + (p i -2Ep + E 2 )p(x,p) 
+ (p 2 -E)$l[V(x)*p(x,p)} - pd x Z[V(x)*p(x,p)} 
--d 2 Jl[V{x)*p{x,p)] - Z[V(x)*pd x p(x,p)} 



+ $s{V{x)*Q[V{x)*p{x,p)}} + $t{V{x)*$t[V{x)*p(x,p)]} 

1 



+ »|v(x)* 



P 



E--d z x \p{x,p) 







(79) 



for x < 0. Clearly, this equation reduces to the one found above when V = 0. 

As a simple check of this result, consider the example of the simple harmonic 
potential V(x) — x 2 . Only the odd-parity wave functions have tp(0) — 0, and 
so only they survive the presence of the infinite potential wall at x = 0. The 
ground state wave function is therefore 

tp(x) = 6(-x) xe-* 2 , (80) 

with energy E — 3. The corresponding Wigner function is p[ip] = 6{—x)p with 

p{x,p) = x 2 erf (x — ip) 7re~ p ~ x 



1 



erf (x + ip)n e p x + x 2 erf (x + ip)ne p x 



-2x(x — ip) 



i\fixp { 



-2x(x — ip) 



r, 2 2 1 2 2 

+ erf(x — ip)p ne~ p ~ x - - erf (x + ip)ire~ p ~ x 

+ ^xe- 2x{x+lp ^ - i^pe- 2x(x+tp) 
+ erf (x - ip)p 2 ne- p2 - x2 . (81) 
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Here erf denotes the error function 

r-X 

erf(x) = / e"* 2 dt . (82) 
Jo 

We have verified that (|81|l does indeed satisfy the equation (|79|) . 

5 Conclusion 

If deformation quantization is truly an autonomous formulation of quantum me- 
chanics, problems that are solved simply using other methods must be treatable 
in it. First-principle deformation quantization of the systems should be possible, 
even if difficult. For stationary states, the Wigner functions should be derivable 
by solving the *-genvalue equation for the Hamiltonian. 

Here we have studied a single quantum particle travelling freely on the x- 
axis, except for the presence of one or two infinite walls. The simple Schrodinger 
wave functions of these systems are well known. The Wigner functions expected 
from the stationary-state wave functions can therefore be easily worked out (see 
p[ip] in eqn {T7J|). Dias and Prata |2J point out, however, that they do not satisfy 
the corresponding *-genvalue equations. 

To cure this problem, they propose the addition of "boundary potentials" to 
the Hamiltonian. They further show that these additional terms are consistent 
with the Schrodinger quantization of the systems. They do not, however, derive 
the extra potentials from first principles. 

In the hopes of filling this gap, we treated the infinite wall and infinite square 
well potentials as limits of solvable, exponential ones. The delta-function well 
is also considered, since it can be treated in a very similar way. Our results do 
not relate easily to the proposal of 0], however. 

A well-defined limit was found for the corresponding *-genvalue equations. 
Happily, this new differential equation is common to all the systems. Further- 
more, the expected Wigner functions p[i/j] satisfy it away from the infinite walls. 
The revised hope is therefore that the new equation can be solved to find the 
correct Wigner functions directly, and boundary conditions on the Wigner func- 
tions can be imposed in the way familiar from Schrodinger quantization. 

That hope has not been realized here. We have not (yet) shown that the 
new differential equations can be solved to find the correct, physical Wigner 
functions. We only verified that the p[?/>]'s satisfied the new equations, outside 
the V = oo regions. 
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Nevertheless, finding a well-defined limit of the *-genvalue equations for 
these systems is a significant first step. The new equation was also generalized 
to allow for a regular potential in addition to infinite walls/wells. 

It might be interesting to use other solvable potentials instead of the expo- 
nential ones. It would at least (hopefully) verify that our results are independent 
of our particular choice of representation of the infinite potential wall as a limit. 
For example, lirn/v^oo x yields the infinite square well. Aspects of this limit 
have already been studied in 



Appendix: Free particle 

In the reviews we have seen, a pure deformation quantization of a free particle 
moving in one dimension is conspicuously absent. Here we attempt such a 
treatment. If nothing else, this appendix should indicate why the free example 
is typically omitted. 

Putting 2m = 1 to get the Hamiltonian H = p 2 , the *-genvalue equation 
becomes ^ 

yflx) P - Ep. (83) 
With h = 1, the imaginary part of this equation is 

pd xP = , (84) 

while 



p'-ll'Df' = E P- (85) 



is the real part. 

The factor p in l|84(l is crucial: we can conclude that when p ^ 0, d x p = 0, 
but not when p — 0. Substituting the ansatz 

P{x,p) = f(p) + S(p)g(x) , (86) 

in (|85|l . we find 

(p 2 - E)f(p) - S( P ) (e + ±%)g(x) = . (87) 



Considering p = yields 

g(x) = 6exp ( 2iVEx ) + b* exp ( - 2iVEx ) , (88) 
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where p = p* has been imposed. Then l|87|l reduces to (p 2 — E)f(p) = 0, solved 

by 

f{p) = a + 5(p-V~E) + a_6{p + ^E), (89) 
with a± arbitrary real constants. 

The terms of (I89|l correspond to plane waves of momentum ±v^E, as can be 
verified by solving 

p*p — p-kp = \J~Ep . (90) 

The expression 

p = a+5{p-VE) + a-5(p+V~E) , (91) 

is the Wigner function of a mixed state of the two momentum eigenstates. The 
terms of H88JI are necessary for coherent superpositions of the two momentum 
eigenstates, and they represent interference between them. That they are re- 
quired can be seen by considering the Wigner function of a simple harmonic 
oscillator in the long-period limit, or the large- width limit of a particle in an 
infinite square well. 

The general result is 




To restrict to pure-state Wigner functions, we impose 

p-kp (x 6(0) p , (93) 

valid for the Wigner functions corresponding to non-normalizable pure states. 
We find 

P*P = 6(0)^(a 2 + + \b\ 2 )S(p- VE) + (a 2 _ + \b\ 2 )5(p + VE) 

+ {a + +a_)5{p) [bexp(2iVEx) + b* exp(-2iVEx)] j . (94) 

The calculation is straightforward, except that we interpret S(p + ^/E)6(p + VE) 
as 5(0)5(p + Ve), e.g., since they yield equivalent results as distributions (i.e. 
upon integration). Similarly, assuming E > 0, we put S(p + ^/E)6(p — \^E) to 
zero, e.g. 

In l|93(l . this result yields the constraint 

\b\ 2 = a+a_ => b = ^Ja+aZe^ , (j) e 3? . (95) 
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The general pure-state solution to the free-particle *-genvalue equation (|83[) is 
therefore 

p = a + 5(p-V~E) + aS(p+y/E) 

+ 2^/o^T S(p) cos (2VEx + </>) . (96) 

On the other hand, calculating 

^ := h /^ e " iP V(^-|)^ + f) , (97) 
with the pure-state wave function 

if) = a+e l ^ x + a_e- %V ^ x (98) 

yields 

p[if>] = \a+\ 2 5{p-VE) + \a_\ 2 S(p + VE) 

+ 5(p)^a* + a-e- 2i ^ x + a+a*_e 2i ^ x ^ . (99) 

Comparing l|99|l with l|96|) reveals a one-to-one correspondence, however, given 
by the relations 

a± = y/a± e#± . - 0+ + 0- = . (100) 

As should be, only the relative phase 0+ — <fi- of (|98|l is relevant to the Wigncr 
function. 
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